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Abstract 



For a power series which converges in some neighborhood of the origin 
in the complex plane, it turns out that the zeros of its partial sums — 
its sections — often behave in a controlled manner, producing intricate 
patterns as they converge and disperse. We open this thesis with an 
overview of some of the major results in the study of this phenomenon 
in the past century, focusing on recent developments which build on the 
theme of asymptotic analysis. Inspired by this work, we derive results 
concerning the asymptotic behavior of the zeros of partial sums of power 
series for entire functions defined by exponential integrals of a certain 
type. Most of the zeros of the n th partial sum travel outwards from 
the origin at a rate comparable to n, so we rescale the variable by n and 
explicitly calculate the limit curves of these normalized zeros. We discover 
that the zeros' asymptotic behavior depends on the order of the critical 
points of the integrand in the aforementioned exponential integral. 

Special cases of the exponential integral functions we study include classes 
of confluent hypergeometric functions and Bessel functions. Prior to this 
thesis, the latter have not been specifically studied in this context. 
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Chapter 1 
Introduction 



The story begins with Gabor Szego, a leading figure in the field of analysis of polyno- 
mials. In 1924 he published a paper [39] in which he examined the radial and angular 
distribution of the zeros of the partial sums — the sections — of the power series for 
the exponential function exp(z). Let 

n f. 
k=0 

be the n th such section. The polynomial s n (exp; z) has exactly n complex zeros while 
the exponential function has none, so Hurwitz's theorem (see, e.g., [23, p. 4]) tells 
us that the zeros of s n (exp; z) must move farther and farther away from the origin 
as n goes to infinity. However, they move in a controlled manner: according to the 
Enestrom-Kakeya Theorem (see Theorem 2.2 in Section 2.1), the zeros of s n (exp; z) 
all lie in the region \z\ < n. Seeing this, Szego studied the behavior of the zeros of the 
polynomials s n (exp; nz). The zeros of these normalized sections all lie in the closed 
unit disk. 

The primary result of Szego's work was that the zeros of the normalized sections 
s n (exp;nz) have as their limit points the simple closed loop 

D = {z e C : \z\ < 1 and \ze^ z \ = 1} . (1.1) 

This curve, and its analogues for other power series, is usually referred to as the Szego 
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-20 20 40 



Figure 1.1: Zeros of the sections s n (exp; z) (n — 1, 2, ... , 70). Lines have been added 
to indicate the zeros corresponding to n = 20, n = 45, and n = 70. 



curve. 

Szego also studied the angular distribution of the zeros. He first showed that the 
mapping w = ze 1 ~ z takes the curve D to the unit circle in the w-plane in a regular 
manner: argw increases monotonically from to 2n as z traverses D from z — 1 in 
the counterclockwise direction. So, if ^(^1,^2) is the number of zeros of s n (exp; z) 
in the sector 61 < arg z < 62 and if z\ and z 2 are the points of D with arguments 9\ 
and 62, respectively, then 

lim »M = w(^)-^i) ^ 

n-5>oo n 2% 



2 



1 Introduction 



0.6 - 



0.4 - 



0.2 - 



0.0 - 



-0.2 - 



-0.4 - 




-0.6 - . . " 

-0.5 0.0 0.5 1.0 

Figure 1.2: Zeros of the normalized sections s n (exp; nz) (n = 1,2,..., 70) with the 
Szego curve D in equation (1.1). 

Essentially this says that, modulo the weight function w, the zeros of the sections are 
uniformly radially distributed. 

Though in this work we are only concerned with the zeros of sections of power 
series, Szego studied the more general question of the roots of the equation 

s n (exp;nz) = Ae nz , (1.2) 

where < A < 1. For A ^ 0, 1 the roots are no longer restricted to \z\ < 1 and may 
accumulate on any part of the curve \ze 1 ~ z \ = 1. For A = 1 the roots accumulate on 
the "arms" of this curve, i.e. the points on \ze 1 ~ z \ = 1 with \z\ > 1. See Figure 1.3 
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for a view of the unrestricted curve \ze 



l-z\ 




Figure 1.3: Part of the curve \ze 



l-z\ 



1. 



Szego also studied the behavior of the analogous question for the power series for 
sine and cosine. By studying the roots of the equation 

(s n (exp; inz) - \e inz ^j + (s„(exp; —inz) - \e~ mz j = 

he was able to apply to these series what he had discovered about the exponential 
series. In particular he deduced that the set of limit points of the zeros of the sections 
s n (sin; nz) and s„(cos; nz) consists of two rotated copies of the part of the limit curve 
D in the right half-plane. Indeed, if D + = D D {z G C : > 0} then the limit 

curve associated with the sine and cosine series is the set 



iD + U -iD + U {x E 



l/e<x< 1/e} . 



Figure 1.4 illustrates the convergence of the zeros of the normalized sections of the 
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cosine function to this curve. 




Figure 1.4: LEFT: Zeros of the sections s n (cos; z) (n = 2,4, . . . , 90). Lines have been 
added to indicate the zeros corresponding to n = 30, n = 60, and n = 90. RIGHT: 
Zeros of the normalized sections s n (cos; nz) (n = 2, 4, . . . , 90) with their Szego curve. 

Most of the results in Szego's paper were rediscovered by Dieudonne [12] in 1935. 

In 1944, Paul C. Rosenbloom showed in his doctoral thesis [31] (and summarized 
in a separate paper [32]) that the behavior described by Szego is in fact a generic 
property of entire functions of positive finite order (see Section 1.3) with a certain 
asymptotic character. To state this result we will require a small amount of notation. 

Let / be an entire function of order < p < oo with 

oo 

f(z) = ^a k z k 

k=0 
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and let 



k=0 

be its n th section. The basic conclusion Rosenbloom came to is that most of the zeros 
of the sections grow on the order of 



Pn = \a n 



-l/n 



As such, he considers the zeros of the scaled sections s n (f;p n z) (another approach 
for determining the appropriate scale factor is outlined in Section 1.3.1). 

Let {iV} be a subsequence of the indices {n} such that the sequence of sections 
{sN{f',z)} has a positive fraction of zeros in any sector with vertex at the origin. 
That is, if ^(#1, #2) denotes the number of zeros of the section s^if] z) in the sector 
9 1 < &rgz < 62, then 

liminf «%M >0 

for any fixed 6\ and 9 2 - Such a subsequence {N} is guaranteed to exist by Theorem 
2.1. Rosenbloom's main result is as follows. 

Theorem 1.1 (Rosenbloom). Suppose that the following conditions hold: 

(1) For some sequence of determinations, f^p^z) 1 ^ converges uniformly to a single- 
valued analytic function g in some subdomain X of the disk \z\ < e 1 ^; 

(2) w = g(z) / z maps X univalently onto a domain X\; 

(3) No limit function of the sequence 

f(pNz) - s N (f;p N z) 



T N z) 



z N 



is identically zero in X ; 

(4) Ttf(z) 7^ in X for N large enough. 

Then the only limit points of the zeros of Sjv(/; Pnz) in X are the points on the 
curve \g(z)/z\ = 1, and their images in Xi under the mapping w = g{z)jz are 
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equidistributed about the unit circle \w\ = 1; that is, the number which accumulate 
about any arc of length a contained in X is asymptotically Na/27r. 

One interesting aspect of Rosenbloom's result is that it allows for the sequence 
fipNz) 1 ^ to converge to different limit functions g in different subregions of the disk 
\z\ < e l l p . One striking example of this is the behavior of the zeros of sections of 
power series for exponential sums of the form 

k=l 

where Cj,Xj G C, as studied by Pavel Bleher and Robert Mallison, Jr. [2]. The 
zeros behave differently in different sectors which are determined by the geometric 
properties of the parameters Xj, as can be seen in Figure 1.5 for the function 

f(z) = 3e (8+2l)z + (-9 + I2i)e^ +7i)z + (2 + z)e ( - 7+4i)z - 5e ( - 6 - 6l)z . (1.3) 




-i.o 



-10 -S 5 10 -1.0 -0.5 0.0 ' 0.5 . 1.0 

Figure 1.5: LEFT: Zeros of the first 100 sections of the exponential sum in equation 
(1.3). RIGHT: Zeros of the first 100 normalized sections with their Szego curve. 

We discuss how the results of Rosenbloom, Bleher, and Mallison relate to the ones 
obtained in this thesis in Chapter 4. 
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1.1 A Nudge Toward Asymptotic Analysis 

James D. Buckholtz was the first to talk about the radial position of the zeros of 
s n (exp;z) in relation to the Szego curve D in equation (1.1). He published a short 
paper [8] on the subject in 1966. His first observation was that all of the zeros of the 
sections lie outside the Szego curve. The proof is so short and elegant that we will 
include it here. 

Theorem 1.2 (Buckholtz). For every positive integer n, neither the curve D nor the 
region it encloses contains a zero of s n (exp; nz). 

Proof. Let zeC with \z\ < 1 and \ze 1 ~ z \ < 1. Then 



1 1 — e nz s n (exp; nz) 



OO U If 

n z 



k=n+l 



k\ 



k=n+l 



oo u 

n 



k\ 



k=n+l 

= 1 — e~ n s n (exp; n) 
< 1. 



Having s n (exp;nz) = here would contradict this inequality. □ 
Buckholtz's result is readily seen in Figure 1.2. 

The proof of this theorem can give a slightly more general result. For n a positive 
integer, if / is analytic in \z\ < n with 

oo 

f( z ) = J2 a kZ k , 

k=0 

let 

n 

s n (f; z) = ^a k z k 

k=0 
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be its n th section. The method of Buckholtz shows that, if s n (f;n)/f(n) > 0, then 
s n (f; nz) has no zeros in the region 



Note that when f(z) = e z we have S n = {z G C : \z\ < 1 and \ze 1 ~ z \ < 1} for all 
positive integers n. 

In addition to describing the direction from which the zeros approach the limit 
curve, Buckholtz used a result from a previous paper of his [7] to examine the rate at 
which they do so. Refining this result would become the central focus of later work 
on the topic. 

Theorem 1.3 (Buckholtz). For every positive integer n, all zeros of s n (exp; nz) lie 
within a distance oflej^fn of D. 

1.2 The Contribution of Newman and Rivlin 

In 1972, Donald J. Newman and Theodore J. Rivlin published a paper [27] in which 
they aimed to establish a zero-free parabolic region for the sections s n (exp; z) de- 
scribed above. However, there was an error in their proof, and so they did not 
actually achieve this goal until their correction [28] was published in 1976. They 
showed that, if c is any positive number satisfying ce c < tt/2, then there is no zero 
in the region {z = x + iy : y 2 < cx} (for related results see, e.g., [33], [35], [34], and 



Though the first paper may not have served its original purpose, the following 
theorem has become important to the theory we're concerned with. 

For this result we will require the complementary error function erfc, defined by 




[36]). 




where the path of integration begins at z and travels to the right to oo. 
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Theorem 1.4 (Newman and Rivlin). For n > 0, define the functions 




where the path of integration begins at w and travels to the right to oo. Then the 
sequence (h n ) converges uniformly to the function 



H(w) = 




on any compact subset of Q(w) > 0, where erfc is the complementary error function. 



The motivation for this result comes from rewriting s n (exp;z) in a form which 
reveals the nature of the "parabolic" arcs of zeros seen in Figure 1.1. First, repeated 
integration by parts will verify that 

s n (exp; z) = J - — e~ l dt. 

Putting z = n + Wy/n and using the substitution ( = w + t/y/n we have 



t (exp; n + Wy /n) _ y/2^(n/e) n 1 f°° f 1 + _C_Y g -^C d(j 



where the path of integration is the horizontal line from w to the right to oo. Note 
that the integral on the right is the function h n (w) defined above. The conclusion of 
Theorem 1.4 can thus be stated as 



g n (exp;n + Wy/n) ^ 1 / 



4 u ; (i.4) 

e n+w^n 2 \ y/2 J 

as n — y oo when w is restricted to a compact subset of ^s(w) > 0. Now if h n has a 
zero w = u + iv with v ^ 0, then z = x + iy = n+(u + iv)y/n is a zero of s n (exp; z) 
which lies on the parabola 

x = {y/vf + u{y/v). (1.5) 
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But if w = u + iv is any zero of the limit function erfc (w/ \/2) in the upper half-plane, 
Hurwitz's theorem tells us that h n will have a zero near w when n is large enough, 
so that s n (exp; z) will have a zero arbitrarily close to the parabola (1.5). In other 
words, the arcs of zeros seen in Figure 1.1 will tend toward parabolas of the form 
x = (y/v) 2 + u(y/v), where w = u + iv is a zero of erfc(w/\/2). 

This behavior can be seen in Figure 1.6. There, w = u + iv is chosen to be the 
smallest zero of erfc(u>/\/2) in the upper half-plane. The parabola associated with 
this zero approximates the upper half of the innermost arc of zeros of s n (exp;z). 




Figure 1.6: Zeros of s n (exp; z) (n — 1,2, ... , 100) and the parabola (1.5) associated 
with the zero of smallest modulus of erfc(w/y / 2) in the upper half-plane. 



The result in equation (1.4) has inspired some similar results for other functions. 
In [18] it is shown that 



s n (E l/x ;R n (l + w 



l + w 



Ei/\(R r , 



erfc(w) 



as n — > oo uniformly when w is restricted to a compact subset of C, where Ei/\ is 
the Mittag-Lemer function of order A and R n is its associated scale factor, both of 
which are described in Section 1.3. A similar result is proved for the ^-functions of 
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order < A < 1, which are described in Section 1.3.3. In [29] it is proved that 

as n — > oo uniformly when w is restricted to a compact subset of C, where ii*i(l; b; z) 
is a confluent hypergeometric function and s n (\F\, b; z) is its n th section — both of 
which are described in Section 1.6. Lastly in [20] we are given that 

when r and n are large with 5 < r/n < 1 — 5 for some S > 0, where 

r 

Br,n{z) = ^ 



fc=0 



for 1 < r < n. 

These results are used to determine the widths of the zero-free regions — the 
"openings" — exemplified in Figures 1.1, 1.7, and 1.10. 



1.3 The Role of the Order of the Limit Function 

Albert Edrei, Edward B. Saff, and Richard S. Varga studied the effect of the order of 
the limit function on the zeros of the sections. In 1983 they published a monograph 
[18] in which they examined the asymptotic character of the zeros of the sections of 
the Mittag-Leffler functions Ei/\ defined by 

CO I, 



where < A < oo. 

Recall that the order of an entire function / — that is, a function which is analytic 
on the entire complex plane C — is defined to be the infimum of all real numbers I for 
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which 

\f(z)\<exp(\zf) 
holds for | z | large enough. Writing 

oo 

f( z ) = Y,akZ k , 

k=0 

we can calculate the order of / directly with the formula 

fclogfc 
i = hmsup- — -— — pr 
k ^oo log(l/|a fc |) 

(see e.g. [3, p. 9] or [37, p. 326]). 

The order of the function E\/\ is seen to be A. As such, E\/\ is called the Mittag 
Leffler function of order A. 



1.3.1 Outline of the Method 

We will outline here the approach used in the monograph. 

The first step is to determine an appropriate scale factor for the sections. Here 
let ^ 

be the coefficient of z k in the power series, and thus let 

n n h 

s n (E 1/x ; z) = J2 a kZ k = WFZTV 

fc=0 fe=0 VA ^ L ) 

be the n th section of the series. We wish to choose the scale factor R n so that the 
inequality 

anK > a k R k n (1-6) 

holds for all nonnegative integers k. For the function Ei/\, the sequence 
{afc_i / 'afcj^x is strictly increasing, so we just need to choose R n to satisfy the in- 
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equality 

— < Rn < —■ (1-7) 
a n a n+l 

Once we do this, inequality (1.6) will in turn be satisfied. By using Stirling's formula 
for the gamma function, we deduce the approximations 



as n — > oo, so if we choose R n such that 



1 fn\ 1 
log^ = x log(-j+-, 



we will indeed satisfy (1.7) for n large enough. 
Define the functions 

U n {z) - 



a n (R n z) 

Rn 



OO oo 
G n (z) = £ — (Rnz)" = J2 ^> 



k=l 



where 



b k (n) = a -^R k n 



for k > —n. By our choice for i? n and properties of the gamma function we see that 
bk{n) — > 1 as n — > oo for any fixed k. It follows that 

GJz) - 



n\-j -. 

1 — Z 



as n — y oo uniformly on compact subsets of the open unit disk. 
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By our definitions we have 

Q n {z) = U n {z) - G n {z). 

As a consequence of the Enestrom-Kakeya Theorem (see Theorem 2.2 in Section 2.1), 
the function Q n has no zeros outside of the open unit disk. If n is large then G n (z) 
is approximately equal to z/(l — z) for \z\ < 1, so if we can find an asymptotic 
representation of the limit function in question, -Ei/a, we can describe the asymp- 
totic character of the zeros of Q n . These zeros are exactly the zeros of the sections 

(El/X'i Rn z )- 

1.3.2 Szego Curves for the Mittag-Leffler Functions 

The definition of the Szego curve for Ey\ is not as easy to state as the one for the 
exponential function; we will need to define it differently in different sectors. The 
curve is defined as the set of all points z — r(9) e td satisfying 

(i) — jx — ® — J\ '■ r (^) is the unique solution of the equation 

r(6) x cos(X6) - 1 - Alogr(fl) =0 

in the interval e" 1 ^ <r(9) < 1, 

(ii) i < < 2vr - i : r{9) = e -V\ 

Thus the Szego curve for Ey^ consists of a circular part and another part whose 
radial component is defined implicitly in terms of its argument. Note that the order 
of the function, A, plays a major role in the definition of the curve. 

1.3.3 Related Results 

In the same monograph, Edrei, Saff, and Varga also derived results similar to those 
above for ^-functions of order < A < 1. Here a function / is said to be an ^-function 
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Figure 1.7: LEFT: Zeros of s n (E\/\\z) (n = 1,2, ...,70) for E w /g (top) and Ei/ 2 
(bottom). RIGHT: Zeros of the normalized sections s n (E\/\\ R n z) [n — 1,2, ... , 70) 
with their Szego curve for £10/9 (top) and Ei/ 2 (bottom). 



if / is entire, /(0) > 0, and 



/« = /«» ft 



where Xk > for all k with 



< 00. 



E 1 

Following Szego's approach in equation (1.2), Natalya Zheltukhina [45] extended 
the work of Edrei, Saff, and Varga on the Mittag-Leffler functions by studying the 
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asymptotic behavior of the roots of the equation 

s n {E\/\\ R n z) = [iEi/\(R n z) 

for \i G C. Later, Zheltukhina and Iossif Ostrovskii [30] performed a similar analysis 
for the classical Lindelof functions, a subclass of the ^-functions studied by Edrei, 
Saff, and Varga. 

1.4 A Careful Study of the Asymptotics 

Five years after the publication of the monograph Varga returned to the problem for 
the exponential function [1 1] , this time with one of his previous PhD students, Amos 
J. Carpenter, and a collaborator, Jorg Waldvogel. 

Their analysis began with the observation that the zeros away from the point z = 1 
approach the limit curve much more quickly than the others. Recall Buckholtz's result 
in Theorem 1.3, that the zeros approach the limit curve D at a rate of 0(l/y/n). 
The first step in the analysis here is to show that this estimate is the best possible 
one when the whole curve is taken into account. Taking full advantage of Theorem 
1.4, the authors prove the following result by tracing the behavior of a zero which 
approaches the point z — 1. 

Theorem 1.5 (CVW). If {z Kn }l =1 are the zeros of s n (exp; nz) and if ti is the zero 
of the complementary error function erfc closest to the origin in the upper half-plane, 
then 

liminf y/H ■ max{dist(^„, D)} > 9fc(ii) + 3f(tx) « 0.636657. 

n— >oo k 

By throwing out the zeros of s n (exp;nz) near z = 1 we should get a different 
estimate for the rate of approach. Indeed, by defining C$ to be the collection of all 
points within a distance S of z — 1, the authors show that the zeros outside of this 
set approach D much more quickly. 

For flCC, define maxdist(fi, C) = sup zen {dist(z, C)}. 
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Theorem 1.6 (CVW). If {z kjn }% =1 are the zeros of s n (exp;nz) and if 5 is any fixed 
number with < 8 < 1, then 

maxdist({;z fcin }£ =1 \ C s , D) = O 

as n oo. 

But the authors noticed something more. As they approach D, the zeros seem to 
lie on regular curves which themselves shrink down to D. Taking inspiration from 
Szego's original analysis, the authors define the intermediate curves 

D n = G C : \ze 1 ~ z \ n = T n v / 2~7m 

where 

n\ _ 1 1 139 

Tn ~ (n/e) n V^n~ ~ + I2n~ + 288n 2 ~ 51840n 3 + " ' 

as n — > oo. Indeed, the curve D n gives a much closer approximation of the zeros of 
•s n (exp;n,z) than does the Szego curve D, as can be seen in Figure 1.8. 

Theorem 1.7 (CVW). If {z Kn }l =1 are the zeros of s n (exp;nz) and if 5 is any fixed 
number with < 5 < 1, then 

maxdist(K n }£ =1 \ C s , D n ) = 0(l/n 2 ) 

as n —7- oo. 

Varga and Carpenter continued their investigation of the asymptotics for the zeros 
of the sections of the exponential series in a second paper [40] . 

Ten years later, Varga and Carpenter published the first of three papers ([41], 
[42], and [43]) in which they carried out a similar analysis of the asymptotics for 
the zeros of the sections of sine and cosine. Following the work of Szego [39], their 
approach essentially began with writing 

2isinz = e lz — e~ lz , 2 cos z = e lz + e~ lz 




1 - z 



< 1, and | arg z\ > cos 



-i 



n 



n 
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-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1.8: Zeros of Si7(exp; 17z) with the intermediate curve Dn and the Szego 
curve D (dotted). 

and 

2is n (sin; z) = s n (exp; iz) — s„(exp; —iz), 2s n (cos; z) = s„(exp; iz) + s n (exp; —iz) 

then applying methods similar to those they used in their analysis of the exponential 
series, though in considerably more detail. 

1.5 A Divergent Power Series 

We shift our focus now to an example of an entirely different sort. In 1996, Karl 
Dilcher and Lee A. Rubel [14] studied the power series 

oo 

J2k\z k . (1.8) 

k=0 

Up to this point the series we have covered have converged on the entire complex 
plane. At the opposite end of the spectrum is this series, which converges nowhere 
but at the origin. It is then incredibly surprising that the zeros of the sections of this 
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series can be wrangled with using essentially the same ideas. Here, the zeros of the 
partial sums 

n 

p n (z) = ^k\z k 

k=0 

move not out to infinity but in toward the origin. So, instead of choosing a scale 
factor which increases with n, one is chosen which decreases with n. 

Theorem 1.8 (Dilcher and Rubel). If a ^ 1 is a complex number, then for all 
positive integers n satisfying 

J 16 4 h ,4 

n > max< — e 1 — a 
- \ 9 

the roots of the equation p n (ez/n) = a lie in the annulus 

1 - r --=<\z\<l + 0-+ a a/—. 

|1 — a\ \Jrt V 7rn 

As a special case of this result, we may use the sharpened form of the Enestrom- 
Kakeya theorem in [1] to see that the zeros of the normalized n th section p n (ez/n) 
all lie in the annulus 

3 

1 = < \z < 1. 

v n 

for n > 97. 




1.6 Inspiration for Studying Exponential Integrals 

The primary inspiration for this current work is a paper by Timothy S. Norfolk [29] 
published in 1998. Studying the confluent hypergeometric functions 



1 F 1 (l;b;z) = T(b)J2 



z k 



where b ^ 1,0, —1, —2, . . ., Norfolk derived results analogous to those of Carpenter, 
Varga, and Waldvogel outlined in Section 1.4. His primary tool was the integral 
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Figure 1.9: LEFT: Zeros of p n (z) = Y%=o k\z k (n = 1, 2, . . . , 70). RIGHT: Zeros of 
the normalized sections p n (ez/n) (n = 1,2, ... , 70) plotted with their Szego curve, 
the unit circle. 



representation 



valid for b > 1. 
If 



1 F 1 (l;b;z) = (b-l) f {1 - t) b ~ 2 e zt dt, 
Jo 



(1.9) 



s n ( 1 F 1 ,6;z) = r(6)5] 



fc=0 



r(A; + 6) 



is the n th section of iFi(l; b; z), we have 

s n (iF 1 ,6;2;) = (6-l) / (1 - t) b - 2 s n (exp; dt, 
where s n (exp; z) is the n th section of the exponential function. Replacing z with nz 
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and subtracting this from equation (1.9) we get the expression 

1 F 1 (l)b',nz)-s n { 1 F 1 ,b',nz) = (b - 1) J (1 - t) fe ~ 2 (e nzt - s n (exp; nzt^jdt 

= {b - 1) / (1 - t) b - 2 e nzt g n {zt) dt, (1.10) 
Jo 

where 

g n (z) = 1 - e~ nz s n (exp;nz). 

This quantity g n (z) was studied by Szego [39], who derived asymptotic approxima- 
tions as n — > oo for z in different regions of the complex plane (see also [11] and [5]). 
In particular he showed that 

where e n (z) = 0(1 /n) as n — > oo uniformly when z is restricted to a compact subset 
of $l(z) < 1. Norfolk applies this approximation of Szego's to obtain an asymptotic 
estimate for the tail (1.10) in the region \z\ < R < 1. 

Norfolk's main result is that the zeros of the normalized sections s n (iFi, b; nz) 
have as their limit points the set 

{z e C : &(z) > 0, \z\ < 1, and \ze l - z \ = l) 
U{zeC: $t(z) = and \z\ < 1/e} 
U {z E C : 9£(z) < and \z\ = 1/e}. 

This curve and some of the zeros can be seen in Figure 1.10. 

We follow essentially the same process in Chapter 3, generalizing Norfolk's result 
to a much larger class of functions. These functions are defined in Chapter 2 and a 
few special cases are listed in Section 4.1. 
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JO 




Figure 1.10: LEFT: Zeros of the sections s„(iF\, 3; z) (n = 1,2,..., 70). RIGHT: 
Zeros of the normalized sections s n (iFi, 3; nz) (n = 1,2, ...,70) with their Szego 
curve. 



Given a power series with finite radius of convergence, Jentzsch's Theorem states 
that every point on the circle of convergence will be a limit point of the zeros of the 
sections of the series. In this sense, the result of Jentzsch is a bridge between the 
results of Dilcher and Rubel in section 1.5 and the results described earlier in the 
chapter for entire functions. We consider a couple simple examples of power series 
with finite radius of convergence in Sections 2.1.1 and 4.2. 

The problem of asymptotic zero distribution has also been treated for a number 
of families of polynomials which are not sections of some power series. We reference 
a few of them here. 

The Daubechies polynomials 



were studied by Djalil Kateb and Pierre Gilles Lemarie-Rieusset [21] and indepen- 
dently by Jianhong Shen and Gilbert Strang [38]. Later, Norfolk and Svante Janson 



1.7 Other Families of Polynomials 
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[20] studied the related polynomials 



k=0 v 7 



where 1 < r < n. The partition polynomials, defined by 



n 



k=l 



where Pk{n) is the number of partitions of n with exactly k parts, were studied by 
Robert P. Boyer and William M. Y. Goh [6]. 
If, in the hypergeometric functions 



we let ai = —n be a negative integer, the series in question are finite and we obtain 
sequences of polynomials. Hypergeometric polynomials of this type been studied 
by various authors (see, e.g., [4], [17], [16], [15], and [46]). Some of these results 
have been used to obtain similar results for the zeros of classical Jacobi orthogonal 
polynomials. 

The problem has also been studied directly for a number of other families of 
orthogonal polynomials (see, e.g., [13], [10], [24], [25], and [5]). Of note are the 
papers of Arno B. J. Kuijlaars, Kenneth D. T.-R. McLaughlin, and Peter D. Miller, 
who apply the methods of Riemann-Hilbert analysis to achieve their results. In 2008, 
Kuijlaars, McLaughlin, Miller, and Thomas Kriecherbauer turned these tools back 
on the original problem of studying the zeros of sections of the exponential series 




and 




where 



T(a + k) 

r(a) 
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[22]. The authors derive full asymptotic series for the zeros in terms of roots of 
appropriately chosen auxiliary equations. 

Though it is only tangentially related to the current theory, we also refer the reader 
to a paper of Christopher P. Hughes and Ashkan Nickeghbali [19] for an interesting 
result about the clustering behavior of zeros of random polynomials. 
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Chapter 2 
Preliminaries 



In this chapter we will discuss the terms and ideas which are relevant to the work. 
Some of these have already been introduced in Chapter 1 and are collected here for 
convenience. 

Our focus will be on the sequence of polynomials given by the partial sums — the 
sections — of a convergent power series. In particular, if Q is some open subset of 
the complex plane C which contains the origin and / : Q — > C is a function which is 
analytic at the origin, then / can be represented by a power series, 

oo 

/(z) = X>z*, (2.1) 

k=0 

in some neighborhood of the origin. We define the n th section of this power series for 
/ to be the polynomial 

n 

s n (f-z) = Y, a kZ k . (2.2) 

k=0 

In later sections we will rely heavily on properties of the exponential function 

OO 

exp(z) =J2h 
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and its sections 

n 

s n {exp;z) 

k=0 

2.1 Strategy and General Tools 

Here we will give a brief description of the method we will use in Chapter 3 to obtain 
our results. We will also collect a few tools for later use. 

Let / and s n (f;z) be as in equations (2.1) and (2.2). For z in the radius of 
convergence of the power series for /, we have 

f{z)-s n {f;z)=t n (f;z), (2.3) 

where 

oo 

t n (f;z) = ^ a k zk 

k=n+l 

is the tail of the power series. Suppose we can find a real, positive sequence {R n } 
such that some subset of the zeros of s n (f; z) grow on the order of R n and let {R n z n } 
be such a sequence of zeros which eventually lies within the radius of convergence of 
the power series. By equation (2.3) this gives us 

/ \Rn^n) ^n(/ Rn^n) 

for n large enough. Theoretically, analytic relations for the scaled zeros z n can 
be derived by determining the asymptotic character of the sequences f(R n z n ) and 
t n (f; R n z n ) as n -> oo. 

An example of this method is given in Section 2.1.1. 

We will need to bound the modulus of the zeros of the sections to simplify the 
estimation of the tails of the power series in Chapter 3. To do this we will invoke a 
particular result of Rosenbloom ([31] and [32]) which will give us precisely the bound 
we need. 

Recall from Chapter 1 that a sequence of sections {sjsr(f\z)} is said to have a 
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positive fraction of zeros in any sector with vertex at the origin if 

liminf i^iA) >Q 

N->oo N 

for any fixed Q\ and 82, where §i?(8i, 82) is the number of zeros of the section s^if] z) 
in the sector 8\ < axgz < 82- Rosenbloom's result is as follows. 

Theorem 2.1 (Rosenbloom). Let 

00 

f(z) = ^a k z k 

be an entire function of finite positive order p and let s n (f; z) be as in equation (2.2). 
Define §° n {R) to be the number of zeros of s n (f;z) in the disk \z\ < R and define 
Pn = |fln|~ 1//n - For any increasing sequence of indices {N} such that the sequence of 
sections {sjv(/; z)} have a positive fraction of zeros in any sector with vertex at the 
origin, the existence of which is guaranteed, and for < r < 1 and e > we have 

rN((e 1/p + e)p N )-r N (pNr) 

liminf — ^ '- > 1 - r p > 0. 

v->oo TV 

Further, the number of zeros of S]y(f]z) satisfying 

\p N z\ > e 1/p + e 



is bounded. 

Part of this result was inspired by the work of Carlson [9] on the angular distri- 
bution of the zeros. 

Many of the recent approaches to the problem (e.g. [29]) have made use of the 
Enestrom-Kakeya Theorem to find this bound. Marden's book [23, ch. 7] is a good 
resource for more information on this theorem as well as for more general results in 
in the same vein. 
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Theorem 2.2 (Enestrom-Kakeya). All zeros of the polynomial 

p(z) = a + a\Z + • • • + a n z n 
having real, positive coefficients aj lie in the ring a < \z\ < (3, where 
a = mm{ak/a k+ i} and (3 = max{afc/a fe+1 } 



for k = 0, 1, . . . , n — 1. 

In [1], Anderson, Saff, and Varga give necessary and sufficient conditions for when 
zeros of p lie on the circles \z\ — a or \z\ — /3. Of particular interest is the maximum 
modulus of the zeros, and as a corollary to their main theorems the authors formulate 
a helpful sufficient condition that the polynomial have no roots on the outer circle 

\z\ = (3. 

Theorem 2.3 (Anderson, Saff, Varga). For p(z) and (3 as defined in Theorem 2.2, 
if f3a\ — ao > 0, then all zeros of p satisfy \z\ < [3. 

The requirement in the above two theorems that all coefficients be positive is 
fairly restrictive. As we mentioned, there are more flexible versions of the Enestrom- 
Kakeya theorem that do not assume this, but they are often not powerful enough 
or simply too difficult to apply to this problem in full generality. For example, the 
best bound these more general theorems can give for the zeros of sections of the 
exponential integrals we will study is precisely twice the bound we need. For this 
reason we instead rely on Rosenbloom's Theorem 2.1. 



2.1.1 A Class of Linear Fractional Transformations 

To illustrate the general strategy described above let us pause to treat a toy example. 

Let ao, A, and B all be positive and real. Define ai = Aa + B and a n = A n ~ 1 ai 
for n > 2. Then the series a kZ k converges to the linear fractional transformation 
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for all z G C satisfying \z\ < 1/A. The sections of this series can be explicitly 
expressed as 

1 - A n z n 

Sn{j;z) = ao + a x z - _ ^ . 

By Theorem 2.3, all zeros of s n (f; z) satisfy \z\ < 1/A. We thus consider the normal- 
ized sections s n (f; z/A), all of whose zeros lie strictly inside the unit circle. 

Given a power series with finite radius of convergence, Jentzsch's Theorem states 
that every point on the circle of convergence will be a limit point of the zeros of the 
sections of the series. Thus every point on the unit circle is a limit point of the zeros 
of the sections s n (f; z/A). If f(z/A) has a zero satisfying \z\ < 1 (which would occur 
at z = —Aao/B), Hurwitz's Theorem (see, e.g., [23, p. 4]) tells us that z = —Aao/B 
will also be a limit point of the zeros of the sections. However, we are only concerned 
with the sections' spurious zeros — that is, sequences of zeros which do not converge to 
zeros of the limit function — so we will exclude these limit points in our calculations. 

We prove the following result. 

Theorem 2.4. Let C be the unit circle, Rs the open region which consists of all 
points within a distance 5 > of the negative real axis, and {zk, n }k=i the zeros of the 
section s n (f] z/A). For flCC, define maxdist(fi, C) = sup^ gQ {dist(z, C)}. Then 

maxdist(K„}£ =1 \^; C) = 0(l/n). 

Proof. Because f(z/A) has no zeros in the set A = {z G C : \z\ < 1} \ Rs, we may 
conclude that the points {z G C : \z\ — l}\Rs are the only limit points in the closure 
of A. 

For z G A we calculate 

s n (f;z/A) = _ ^ n+ i ai 
fiz/A) Aa + Bz' 

So, if z is a zero of s n (f] z/A), it must satisfy the relationship 

z"* 1 = (Aa + Bz)/ ai . (2.5) 
For z G A it is clear that there exist positive constants C% and Ci such that 
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G\ < \(Aao + Bz)/ai\ < C2, so that on taking absolute values and (n + l) th roots in 
(2.5) we find that the spurious zeros of s n (f; z/A) satisfy 

\z\ = l + 0(l/n), 

as desired. □ 
This result is illustrated in Figures 2.1 and 2.2. 




Figure 2.1: Zeros of s n (f; z) (n = 1, 2, . . . , 60), where / is the linear fractional trans- 
formation in (2.4) with do = A = B = 1, and their Szego curve, the unit circle. 

2.2 The Exponential Integral Functions 

The functions we will study will be defined by integrals of the form J^ a (p(t)e zt dt. 
The restrictions we will place on the function (p are determined essentially by the 
abilities of Watson's Lemma, which are discussed in the next section. 
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1.0 - 




Figure 2.2: Zeros of %(/; z), where / is as in (2.4) with a — A — B — 1, with the 
intermediate curve (the modulus of equation (2.5)) and the unit circle. 

Suppose < a, b < oo and let <p : [—a, b] — > C U {oo} be a function satisfying 




\<p(t)\dt < oo 



and <p(t) = (t + aYf^t + a) = (b - t) u f 2 (b - t), where 

(1) fi, v e C with 9fc(/z) > -1 and > -1, 

(2) fx, f 2 : [0, a + b] — > C U {oo} with /i(0) and /2(0) both finite and nonzero, 

(3) in a neighborhood of t = 0, both f[(t) and f' 2 (t) exist and are bounded. 
Define 

= / ip(t)e zt dt. 
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The function F is entire, and its sections are given by the formula 




Notice that F(z) is the exponential generating function of these particular integral 
moments of (p. 

Properties (1) and (2) above essentially serve to ensure that the integral F(z) 
converges. Property (3) allows us to determine a simple error term in the asymptotic 
expansion of F in Watson's Lemma. 

An important example of functions of this type are the Bessel functions of the 
first kind J a with 3fi(a) > —1/2. Indeed, in this case we have 



We discuss in detail how the main result of Chapter 3 applies to these Bessel functions 
in Section 4.1. 

2.3 Watson's Lemma 

The primary tool in this work is Watson's Lemma. For a thorough discussion of this 
result in a general setting see [26]. 

In the following, A is a complex parameter. 

Theorem 2.5 (Watson's Lemma). Suppose < T < oo and <p : [0, T] — > C U {oo} 
is a function satisfying 



Jo 

and ip(t) = t a h(t), where dt(a) > —1, h(0) ^ 0, and h'(t) exists and is bounded in a 
neighborhood oft — 0. Then the exponential integral 
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is finite for all 5R(A) > ; and 



$(A) 



fe(0)r(cr + l) 



+ 0(A- CT - 2 ) 



as A — > oo | arg A| < for any fixed < 9 < n/2. 

Though this form of Watson's Lemma only gives an asymptotic for $(A) as A — > oo 
to the right, it can easily be extended to address the case when A — > oo to the left. 
However, to do so we must assume that T is finite. 

Corollary 2.6. Suppose < T < oo andip : [0,T] — > CU{oo} is a function satisfying 



and ip(t) = (T - t) a h(T - t), where 3?(<r) > -I, h(0) ^ 0, and h'(t) exists and is 
bounded in a neighborhood oft — 0. Then the exponential integral 





is finite for all 3?(A) > ; and 



$(A) 



/gV+i) 



e TX + 0(\-°- 2 e TX ) 



as A — > oo with | arg A| < for any fixed < 6 < tt/2. 



Proof. We have 




fe(0)r(g + l) 
A-+ 1 



+ 0(A— 2 ) 
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as A — > oo with | arg A| < 9 for any fixed < 9 < tt/2, by Watson's Lemma. □ 

Watson's Lemma and its corollary will be used in Chapter 3 to determine the 
asymptotic character of the exponential integrals described in Section 2.2 as well as 
the coefficients of their power series. 
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Chapter 3 
Main Results 



Recall from Section 2.2 that we are concerned with functions of the form 



F(z) 



ip(t)e zt dt 



with ip satisfying some light requirements. The n section of F is the polynomial 



Z f 

k=o J - a 



^{t)t k dt. 



The statement and proof of the main result depends on the relative sizes of a and 
b and of 3?(/x) and To this end, define c = max{a, b} and 



^t(n) if a > b, 

&(y) if a < b, 

min{3fj(/i),3?(z/)} if a = b. 



Lastly, let {iV} be an increasing subsequence of the indices {n} for which the 
sequence of sections {s^(F; z)} has a positive fraction of zeros in any sector with 
vertex at the origin (see Section 2.1). If a = b and dt(fi) = dt(u), we also impose the 
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condition that the indices {N} are chosen so that quantity 

(-1)^(0)^ + 1) + f 2 (0)T{u + lK-W" (3.1) 

is bounded away from 0. Such a subsequence is guaranteed to exist by Theorem 
2.1. The condition in (3.1) ensures that we can use the asymptotic representations 
derived in Lemmas 3.4 and 3.7 without incident. 
The main theorem is as follows. 

Theorem 3.1. It is true that 

(i) Every point on the curve 

D atb = jz G C : &(z) < 0, \z\ < ^, and \cze 1+az \ = 1 

' - 1 l-bz\ 



U iz G C : > 0, \z\ < -, and \cze L - bz \ = 1 

is a limit point of the zeros of the sections sn(F; Nz). 

(ii) The only limit points of the zeros of the sections s^(F; Nz) on the imaginary 
axis are those on the line segment 

Dimag = |-2 G C : = and \z\ < — 

The limit points of the zeros of F(Nz), if there are any, are a subset of the 
imaginary axis. If every point on the imaginary axis is a limit point of these 
zeros, then every point on D imag is a limit point of the zeros of the sections 
s N (F;Nz). 

(Hi) Let {zn} be a sequence of complex numbers such that sn{F; Nzn) = for all 
N and such that the sequence has a limit point in the region 

z G C : < and z ^ -- 

a 
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Then the elements of the sequence satisfy 

\cz N e 1+az "\ = 1 + (t - Sty,) + iy^L + 0(1/N) 

as N — > oo. 

(iv) Let {zn} be a sequence of complex numbers such that sn(F; Nzn) = for all 
N and such that the sequence has a limit point in the region 

1 2 e C : > and z^^\. 

Then the elements of the sequence satisfy 

\czn*-*"\ = 1 + (t - »w + + 

as N oo. 

It is interesting to note that the zeros will eventually approach the limit curve 
D a iy from the inside or from the outside depending on the signs of the quantities 
f - &(//) + 1/2 and f - + 1/2. For example, if ^ - 3?(//) + 1/2 > then the zeros 
will eventually approach the part of D a)0 in the left half-plane from outside the curve. 
If either of these quantities is zero then this theorem does not give any information 
about the direction from which the zeros approach the relevant part of the curve. 

Figures 3.1 through 3.6 showcase the limit curve D a) b U -D imag and the zeros of the 
sections for a few exponential integrals. 

The remainder of this chapter is dedicated to the proof of this theorem. We begin 
by finding two asymptotic estimates we will require in our calculations. 

Lemma 3.2. As n — > oo, 

F(nz) = A(0)r(/i + l)(-7iz)-^ 1 e- an * (l + 0(l/n)) 
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when z is restricted to a compact subset of 9fi(z) < ; and 

F{nz) = f 2 (<d)T(v + l)(nz)- v - l e bnz (l + 0(l/n) 
when z is restricted to a compact subset of dt(z) > 0. 

Proof. This follows from a direct application of Corollary 2.6. To see this, suppose 
first that z is restricted to a compact subset of < 0, and make the substitution 
t = b — s in the integral for F(nz) to get 

/b pa+b 
<p(t)e nzt dt = e bnz / <p(b - s)e~ nzs ds, 
■a JO 

which, after replacing z with —z, is of the form required by the corollary. Next, sup- 
pose that z is restricted to a compact subset of > 0, and make the substitution 
t = s — a in the definition of F(nz) to get 

/b pa+b 
<p(t)e nzt dt = e~ anz / cp(s - a)e nzs ds, 
■a JO 

which is also of the required form. □ 

Lemma 3.3. We have 
rb 

<p(t)t n dt = (-l) n A(0)r(/i+ l)n-"- 1 a n+ " +1 + 0(n-^~ 2 a n ) 
+ / 2 (0)r(i/ + l)n- u - 1 b n+u+1 + 0{n- u ~ 2 b n ) . 

as n oo. 

Proof. If a 7^ we calculate 

pO 

ip(t)t n dt = (-a) n / ip(t)e nlog{ - t/a) dt 



{-a) n / (p(s - a)e nlog(1 " s/a) ds 
Jo 

{-a) n [ s^/!(s)e nlog(1 - s/a) rfs. 



o 
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Letting s = a(l — e r ) gives 

/O poo 
ip(t)t n dt = (-l) n a" +/i+1 / (1 - e- r Yh{a - ae- r ) e - r e~ nr dr 
a JO 

POD 

= (-l) n a n+ ^ +1 / r^ a (r)e~ nr dr, 
Jo 

where 

Mr)=( h(a - ae- r )e- r 



r 

has a bounded derivative in a neighborhood of r = 0. We may now apply Watson's 
Lemma to conclude that 

/ tp(t)t n dt = (-l) n V' a (0)r( At + l)n- /1 - 1 a n+ ' i+1 + O (n~^V) 

J — a 

= (-l) n /i(0)r(// + l)^-^ 1 ^^ 1 + 0(n-^ 2 a n ) . 
Using an identical argument we find that 

/ ip(t)t n dt = f 2 {Q)T{v + l)n- v - l b n+v+l + 0(n- v - 2 b n ) , 
Jo 

which completes the proof. □ 

The actual asymptotic character of the integral in the above lemma depends on 
the relative sizes of a and b and of 3?(/i) and so for convenience we'll organize 

the possible outcomes in a separate lemma. 

Lemma 3.4. If a > b then 

f <p(t)t n dt = (-l) n /i(0)r(^ + 1) n^V^ 1 (l + 0(l/n)) , 

J — a 

if a < b then 

J" V {t)t n dt = / 2 (0)r(i/ + 1) n- / - 1 b n+v+1 (l + 0(l/n)) , 
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and if a = b we have three cases: 

(i) ifR(p) > &(u) then 

f <p(t)t n dt = f 2 (0)T(v + I) n-'-'a^ 1 (l + O^^-^) +0(l/n)\ 

J —a 

(ii) ifft(n) < then 

j <p(t)t n dt 

J —a 

= (-l) n+1 / 1 (0)r( / u + l)n^- 1 a n+ ^ +1 (l + O(n^ ) - 5RM ) +0(l/n)), 
(Hi) and if 3ff(/x) = dt(u) then 

j ip(t)t n dt = Gi(a, fx, u, n, z) n-»- l a n+ » +1 (l + 0(l/n)) , 

./ —a 

Gx(a, /i, !/, n) = (-l)" +1 /i(0)r(/i + 1) + / 2 (0)r(i/ + 1)0"-^", 

ifGi(a,(j,, u,n) ^ 0, 
as n — t- 00. 

Proof. This follows directly from the statement of Lemma 3.3. □ 

3.1 Restricting the Zeros 

In this section we will prove the following lemma. 

Lemma 3.5. The limit points of the zeros of the sections sn{F;Nz) lie in the disk 
\ z \ < l/ c ? where c = max{a, b}. 

From Stirling's formula 

n! ~ ^— j V2nn 
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3.2 Gathering the Formulas 



we have 

(nl) 1 ^ ~ - 

e 

as n — > oo, and with the aid of Lemma 3.4 we calculate 



ip{t)t N dt 



-l/N 



C 



as N — > oo, where the subsequence of indices {iV} is as defined above Theorem 3.1. 
Combining these we see that the power series coefficients of F(z), which are given by 



1 



an = ~. <p(t)t n dt, 
n\ 

satisfy 

I i-i/iv N 

Pn = \a N \ 1 , 

ec 

where is as defined in Theorem 2.1. Thus the order p of F is calculated to be 

logn 

p = hmsup = 1, 

n^oo log p n 

where all indices n are taken into account (see, e.g., [3, p. 9] or [37, p. 326]). 
By Theorem 2.1 we have, for any e > 0, 



^((l + e)AT/ C 
liminf = 1, 

N-+oo N 

where ^° N ((1 + e)N/cj is the number of zeros of sn(F; z) in the disk \z\ < (1 + e)N/c. 
From this we conclude that the limit points of the zeros of the sections sn{F; Nz) lie 
in the disk \z\ < 1/c, as desired. 

3.2 Gathering the Formulas 

By definition we have 

b 

nzt 



F{nz) = / ip(t)e nzt dt 
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3.2 Gathering the Formulas 



and 

r-6 



s n (F;nz) = j ip(t)s n (exp;nzt) dt. 
Subtracting these we get 



F(nz) - s n (F; nz) = j <p(t) (e nzt - s„(exp; nzt)) dt 

V (t)e nzt g n (zt)dt, (3.2) 



where 

g n {z) = 1 - e~ nz s n (exp;nz). 
It was shown by Szego in [39] (see also [11], [5], and [29]) that 

where e n (z) = 0(1 /n) as n — > oo uniformly when z is restricted to a compact subset 
of < 1. Upon substituting this into equation (3.2) we get 

F(nz) - s n (F; nz) = ^= f t n+1 (l - e n (zt)) dt. (3.3) 

\/2nn J- a 1 — zt V / 

Our next step is to estimate this integral. The following lemma is proved by an 
argument similar to the one used to prove Lemma 3.3. 

Lemma 3.6. 

f b JfR- t n+1 dt = (_i)"+i /i(°) r fc + V n -M-i a n+M+2 + (n-»- 2 a n ) 
J_ a 1 - zt y ' 1 + az v ; 

l-bz v ' 

as n — > oo uniformly when z is restricted to a compact subset of the doubly-punctured 
plane {z 6 C : z ^ —1/a and z ^ 1/b}. 

As in Lemma 3.3, the actual asymptotic character of the integral in the above 
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3.2 Gathering the Formulas 



lemma depends on the relative sizes of a and b and of 9?(/i) and $l(v), so we will again 
organize the possible outcomes in a separate lemma. 

Lemma 3.7. If a > b then 

J_ a l- zt v 7 1 + V v 7 V 



if a < b then 

fb 



£ J^L ( »« * = /»(Q>r_fr + 1) „-.^ (i + 0( i/„)) , 



and if a = b we have three cases: 
(i) ifU{n) > £/ien 



r+ i ^ = / 2 (o)r(»/ + i) n -,-i a n + , +2 A + / rm-rm) + 0{l/n) \ 

J_ a 1 — zt 1 — az V / 

f^Lt^dt 

J -a 1 " 



= 1) + ( W »M-»M) + Q(l/n)), 

fnij and i/3?(/i) = 9ft(z/) £/ien 

G 2 (a,n, u,n,z) 

= (-l) n+1 /i(0)r(/i + i)(i _ az ) + / 2 (0)I> + 1)(1 + az)a l/ -^-", 
ifG 2 (a,n, v,n,z) ^ 0, 
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3.3 Concluding the Argument 



as n — » oo uniformly when z is restricted to a compact subset of the doubly-punctured 
plane {z 6 C : z / —1/a and z 7^ 1/6}. 

Proof. This follows directly from the statement of Lemma 3.6. □ 
As a consequence of these two lemmas, equation (3.3) becomes 

when z is a zero of the scaled section s n (F;nz) in the region |9ft(.z)| < 1/c, where 
c = max{a, b}. 

3.3 Concluding the Argument 

Let {N} be a sequence of indices as defined above Theorem 3.1. We showed in 
Section 3.1 that the limit points of the zeros of the sections sn{F; Nz) lie in the disk 
\ z \ < 1/c, where c = max{a, b}. Define 

A = < z E C : \z\ < -, z ^ , and 2 ^ - 

L cab 

Suppose first that {zn} is a sequence in C such that Sjv(F; Nzjsi) = for all iV 
and such that the sequence has a limit point in A fl {z G C : !R(z) < 0}. Note that 
this implies there is a 5 > such that \z^ + 1 / a\ > 5 for N large enough. 



It follows from Lemma 3.2 that 



\F(Nz N )\^ N = \e~ aZN \ (1 - ORGu) + 1) ^ + 0(1/^)) (3.5) 



as iV — > 00, and if 

f 

3?(/x) if a > b, 

if a < b, 

min{3?(/x), 3?(i / )} if a = b, 
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3.3 Concluding the Argument 



then we have from Lemma 3.7 that 

ecz N \ (l - ^ + ^ + 0(1/JV)) (3.6) 

as — )• oo. Upon substituting equations (3.5) and (3.6) into equation (3.4) we see 
that these zeros zn satisfy 

\cz N e 1+az -\ = 1 + (£ - Rfo) + ^ + 0(l/iV) 

as iV — > oo, which proves part (iii) of Theorem 3.1. 

Suppose now that {z^} is a sequence such that s^(F; Nz N ) = for all iV and 
such that the sequence has a limit point in A fl {z G C : 3?(z) > 0}. Note that this 
implies there is a 5 > such that \zn — l/b\ > 5 for N large enough. 

Here it follows from Lemma 3.2 that 

\F{Nz N )\^ N = |e te -| (l - (»(»/) + 1) ^ + 0(l/iV)) (3.7) 

as A" — > oo. Subtituting this and equation (3.6) into equation (3.4) we see that these 
zeros zn satisfy 

Ic^e 1 -^! = 1 + (t - »(!/) + ^ + 0(1/N) 

as A" — > oo, which proves part (iv) of Theorem 3.1. 

We have so far shown that the limit points of the zeros of the sections sn(F; Nz) 
with ^ft(z) must lie on the curve D a ^ as defined in part (i) of Theorem 3.1. That 
every point on D a ^ is such a limit point follows from the choice of the subsequence 
{A^} which ensures that the sequence of sections {sn(F; z)} has a positive fraction 
of zeros in any sector with vertex at the origin. It is straightforward to show that for 
any < 9 < 2ir there is a unique r > such that re %e G D a ^. This proves part (i) of 
Theorem 3.1. 

Finally we will prove part (ii). From the asymptotic expansion for F(nz) in 
Lemma 3.2 we see that the limit points of the zeros of F(nz) must lie on the imaginary 



e N z N+1 



<P(t) 

V2nN J - a I- zt 



lN+1 



dt 
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3.4 A Few Examples 



axis. 

If {zn} is a sequence of complex numbers such that sn(F; Nzn) = for all N and 
such that the sequence has a limit point in {z G C : $l(z) = and $S(z) > l/(ec)} 
then by equation (3.4) and Lemma 3.7 we must have F(Nzn) oo. However, if 
j/eM then 



\F(iny)\ 



b 



ip(t)e inyt dt 



b 



< / \<p(t)\dt, 



so that such a sequence of zeros cannot exist. Hence any limit points on the imaginary 
axis must satisfy ^s(z) < l/(ec). 

Further, if {z^} is a sequence of zeros which has a limit point in \z\ < l/(ec), 
then by equation (3.4) and Lemma 3.7 we must have F(Nzn) — > 0. In other words, 
the zeros {z n} must approximate the zeros of F(Nz). Thus part (ii) is proved, 
completing the proof of Theorem 3.1. 

3.4 A Few Examples 

In this section we'll look at what Theorem 3.1 says about a few example exponential 
integrals. 

In Figures 3.1 and 3.2 we plot the zeros of the sections and of the normalized 
sections, respectively, of the function 

r3/2 

F^z) = J (3/2 -t) l l 2+i e zt dt. 

In the notation of Theorem 3.1 we have a = 2, b = 3/2, /i = 0, and v = 1/2 + i. The 
Szego curve associated with the normalized zeros is the set 

^2,3/2 U Amag = j* G C : &(z) < 0, \z\ < i and \2ze 1+2z \ = 1 j 



\z\ < -, and 


2ze 1 ~l z 


-} 









U <J z e C : ft(z) = and \z\ < — 
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3.4 A Few Examples 



Zeros of sjv(-F\; Nz) which converge in ^ft(z) < to a point different from z = —1/2 
satisfy 

\2ze^\=l+ l ^f + 0(l/N) 

and thus eventually approach the curve D 2 ^/2 from the outside. Zeros which converge 
in $l(z) > satisfy 



1 + OQ./N). 



In Figures 3.3 and 3.4 we plot the zeros of the sections and of the normalized 
sections, respectively, of the function 



F 2 (z) = j 1 (1 - t)\t + i)-i/2-2i e *t dt 



In the notation of Theorem 3.1 we have a — b — 1, fi — —1/2 — 2i, and v — 4. The 
Szego curve associated with the normalized zeros is the set 

£>i,i U Amag = {z e C : U(z) < 0, \z\ < 1, and |^e 1+z | = 1} 

U {z E C : SRO) > 0, \z\ < 1, and |ze^ 2 | = 1} 
U {z E C : = and \z\ < 1/e} . 

Zeros of sn{F2] Nz) which converge in 9ft(.z) < to a point different from z — — 1 
satisfy 

= 1 + ^ + 0(1/^) 

and thus eventually approach the curve .D^i from the outside. Zeros which converge 
in > to a point different from z = 1 satisfy 

|«'-| = l-^ + O(l/JN0 

and thus eventually approach from the inside. 

In Figures 3.5 and 3.6 we plot the zeros of the sections and of the normalized 
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sections, respectively, of the function 



/ 19/36 
(t- 1/2) V* eft. 
-17/36 

In the notation of Theorem 3.1 we have a = 17/36, b = 19/36, and fj, = v = 0. The 
Szego curve associated with the normalized zeros is the set 



-^17/36,19/36 U Amag 



36 

z e C : < 0, \z\ < — , and 

19 



— 2;e + 36 z 




36 





U Iz E C : »(z) > 0, \z\ < 



36 
19 : 



U <z e C : SR(^) = and \z\ < 



and 



36 
19e 



19 i_ii 2 




36 





Zeros of s^(F 3 ;Nz) which converge in %t(z) < satisfy 



19 i+i? 2 
— + 36 z 

36 



logiV 
l + ^T + 0(l/iV) 



and thus approach the curve 1^17/36,19/36 from the outside. Zeros which converge in 
3t(z) > to a point different from z = 36/19 satisfy 



19 119. 

— ze 36* 
36 



logiV 

i + ^r + o(i/A0 



and also approach -Di7/ 36j i9/ 36 from the outside. 
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3.4 A Few 



Examples 
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3.4 A Few Examples 
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Chapter 4 



Discussion 



Rosenbloom's main theorem, Theorem 1.1, is quite general. For instance, it can be 
used in conjunction with Theorem 2.1 to determine the limit curve described in part 
(i) of our Theorem 3.1. The new result in this current thesis is the derivation of the 
asymptotic order at which the zeros approach this limit curve as well as information 
about the direction at which they do so in most cases. This information is given in 
parts (hi) and (iv) of Theorem 3.1. 

In light of Theorem 1.1, one consequence of Lemma 3.2 is that the zeros of the 
sections of the exponential integrals we have studied have the same limit curves as 
the zeros of particular exponential sums. Indeed, if a, b > and 

F(z) = I ^(t)e zt dt, 

J —a 

where ip is an appropriate extended-complex-valued function as described in Chapter 
2, the zeros of the sections s n (F; nz) have the same limit curve as the zeros of the 
sections of the sum 

e~ az + e bz . (4.1) 

Zeros of sections of exponential sums such as these were previously studied by Bleher 
and Mallison [2] (see Figure 1.5 and the paragraph immediately following Theorem 
1.1). 

One interesting aspect of the current work is that the asymptotic rate of approach 
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4.1 Special Cases of the Exponential Integrals 



of the zeros of the sections s n (F;nz) depends on the order of the critical points of 
ip(t) at t — —a and t = b (see parts (iii) and (iv) of Theorem 3.1). In most cases 
the zeros approach the limit curve at a rate of clogn/n for some nonzero constant c 
(either positive or negative), and in other cases at a rate of 0(l/n). This stands in 
contrast with the work of Bleher and Mallison, who showed that the zeros of sections 
of exponential sums, such as the one in equation (4.1), always approach the arcs of 
the limit curve at a rate of clogn/n for some positive constant c. 



4.1 Special Cases of the Exponential Integrals 



with b > 1 studied by Norfolk [29] are a special case of the exponential integrals 
studied in this thesis. Our result extends some of Norfolk's results to the case of 
complex b with 9ft (6) > 1. A few of his results, notably the asymptotic rate of 
approach of the zeros to the "corners" of the Szego curve, were not replicated and 
will be treated in future work. 

Another important special case of the exponential integrals is the class of Bessel 
functions of the first kind 



which is of the required form once z is replaced with — iz. As far as we can tell, 
asymptotics for the zeros of sections of the Bessel functions have not been previously 
studied. We state the result formally as a corollary to Theorem 3.1. This corollary 



As discussed in Section 1.6, the confluent hypergeometric functions 





For 9ft (a) > —1/2 we have Poisson's integral representation (see, e.g., [44]) 
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4.1 Special Cases of the Exponential Integrals 



is illustrated in Figure 4.1. 

Corollary 4.1. Let J a be the Bessel function of order a with 9ft(a) > —1/2, and for 
n even let 

z a n ^ 2 ( l) k z s "' 

Sn{J a \z) = {^) E r (Hl)r(Ha + l) (2, 

fre z£s n th section. Then, for a subsequence {N} of the indices {n} as defined in 
Theorem 2.1, the zeros of the normalized sections s^{J a ;Nz) have as their limit 
points the set 

D(J a ) = {z G C : S(z) > 0, |2| < 1, and \ze 1+iz \ = 1} 

u{zGC:3(z)<0, \z\ < 1, and \ze l ~ iz \ = l) 
U {x 6 R : - 1/e < x < 1/e}. 

Conversely, every point of D(J a ) is a limit point of zeros. Zeros of s^{J a \Nz) in 
^s(z) > which do not converge to the point z = i satisfy 

K-|=l + f^ + O(l/A0 

as N — >■ 00, and zeros m ^(z) < which do not converge to z = —i satisfy 

1 ■ , logiV 
\z^-"\ = l + -^- + 0(l/N) 

as N 00. Furthermore, if a is real, the modified indices {N} in this result may be 
replaced everywhere by the original even indices {n}. 

Proof. The only parts of Corollary 4.1 which do not follow immediately from Theorem 
3.1 are the claims 

(i) Every point of the line segment {16R: — l/e<a;<l/e}isa limit point of 
zeros, and 

(ii) If a is real, the modified indices {iV} may be replaced by the original indices 

w. 
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4.1 Special Cases of the Exponential Integrals 



To prove claim (i) we cite the known result [44] that the large zeros of the Bessel 
function J a are given by the asymptotic expansion 

k+ a 1\ 4a 2 -1 (4a 2 - 1) (28a 2 -31) 



2 47 8 (* + f-i)7T 384 [(A;+f-i) 7 r] 2 

Every point of the imaginary axis is thus a limit point of the zeros of J a (Nz), so that 
every point of the segment — l/e<x<l/e}isa limit point of the zeros 

of the normalized sections s^iJa] Nz). Thus claim (i) is proved. 
To prove claim (ii) we write 



s n (J a ;-inz) = Pn{z 2 ) 



where 



n/2 



Pn(z) = J2 



n 



2k 



Z 



k 



, A k T(k + l)T(k + a + 1) 
By the Enestrom-Kakeya theorem (Theorem 2.2), all zeros of P n (z) satisfy 

2n + 4 T(f + a + 2) 



z < 



n 2 



r(f + a + iy 

We apply Stirling's formula for the Gamma function to see that 

r(f + a + 2) 



n 



1 



r(f + a + l) 2' 

whence 

2n + 4 T(f + a + 2) 
n 2 'r(f + a + l) " 

as n — > 00. This tells us that the limit points of the zeros of the polynomials P n , and 
hence of the sections s n (J a ] nz), lie in the closed unit disk. This replaces the use of 
Theorem 2.1 in the proof of Theorem 3.1, which is the origin of the restriction of the 
indices to the subsequences {N}. This proves claim (ii), completing the proof of the 
corollary. □ 
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-40 -20 20 40 -1.0 -0.5 0.5 1.0 



Figure 4.1: LEFT: Zeros of the sections s n ( Jo; z) (n = 2, 4, . . . , 90). RIGHT: Zeros 
of the normalized sections s n ( Jo; nz) (n = 2, 4, . . . , 90) with their Szego curve. 

4.2 Generalizing the Method 

Though most of this thesis has dealt with sections of entire functions, we showed 
in Section 2.1.1 that it is sometimes possible to perform this analysis on sections of 
series with positive finite radius of convergence. As another example, if 

1 oo 

/(*) = 7T^yi =£(* + !)**> ( 4 - 2 ) 

we can use a similar approach to show that the zeros of the sections s n (f; z) which 
do not converge to the point z — 1 satisfy 

loe; n 

\ z \ = i _^i + o(l/n) 
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as n — > oo (see Figure 4.2). Usually this analysis will require the use of the Enestrom- 
Kakeya theorem or one of its variants to ensure that the zeros eventually lie within 
the radius of convergence, but somtimes tricks can be employed if this is not actually 
the case. If, say, the relevant zeros all lie outside the radius of convergence, we could 
try studying the zeros of the polynomials z n s n (f; 1/z) instead. 




Figure 4.2: Zeros of s n (f; z) (n = 1, 2, . . . , 75), where / is the rational function in 
(4.2), and their Szego curve, the unit circle. 

In the future we would like to investigate the possibility of studying exponential 
integrals of integer order greater than one. For example, it would be interesting to 
see if we could simply replace z with z m for some integer m and derive analogous 
results for the new function 

F m {z)=F{z m )= f V (t)e zmt dt, 

J —a 

which would have exponential order m. 

Lastly we speculate that it may be possible to iterate the process described in 
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4.2 Generalizing the Method 



this work. If / is an entire function amenable to the current method — that is, if 
reasonably-detailed asymptotics for f(z) and its tail t n (f;z) = f(z) — s n (f;z) are 
known — then it may be possible to determine useful asymptotics for a new function 
f(z), defined by 

f(z)= I <p(t)f(zt)dt, 

J a 

as well as for its tail t n (f; z). The zeros of the sections of the new function / could 
then be studied. It would be interesting to see what, if any, characteristics of the 
zero distribution associated with / carry over to the one associated with /. 
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Appendix A 

Mathematica Code for the Plots 



Here we give examples of the code used to create the various plots in this thesis. All 
plots were created in Mathematica 7. 

The zeros of the sections of any function f(z) which is analytic at z = can be 
plotted using code like the following. 



f [x_] :=E~x; 




s [ n_, x_] : =Normal [Series [ f [ z ] , 


{z, 0,n}] ] /.z->x; 


numpolys=2 ; 




start=l ; 




allzeros= { } ; 




For [k=start, k<numpolys + start , 


k + +, 


newzeros=x/ .NSolve [ s [k, x] 


==0,x, 70] ; 


For [ j=l, j<=Length [newzero 


s] , 


AppendTo [all zeros, 




{Re [newzeros [ [ j ] ] 

] ; 

] ; 


] , Im [newzeros [ [ j ] ] ] } 


] ; 

Show [Graphics [ {Point [allzeros 


] }] ] 



Code Snippet 1. 
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A Mathematica Code for the Plots 



Here, the function in question is f(z) = e z , and the zeros of the first 20 (from start 
to numpolys+start-l) sections are plotted, producing an image like Figure A.l. 



Figure A.l: Output of Code Snippet 1. 
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A Mathematica Code for the Plots 



To plot the zeros of the normalized sections, simply replace x with k x in the 
NSolve function. The modified block from Code Snippet 1 follows. 

numpolys=2 ; 
start=l ; 
allzeros= { } ; 

For [k=start, k<numpolys + start , k + +, 

newzeros=x/ .NSolve [s [k, k x] ==0, x, 70] ; 
For [ j=l, j<=Length [newzeros] , 
AppendTo [allzeros, 

{Re [newzeros [ [ j ] ] ] , Im [newzeros [ [ j ] ] ] } 

] ; 

] ; 

] ; 

Code Snippet 2. 







•••• 



Figure A. 2: Output of Code Snippet 2. 
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A Mathematica Code for the Plots 



The Szego curves can be drawn with the ContourPlot function. The code to 
draw the Szego curve for the exponential function is given below. 



Show [ 






Graphics [ {Point [al 


lze 


ros] } ] , 


ContourPlot [ 






Abs [ (x+I y) E" 


(1- 


(x+I y) ) ]==1, 


{x,-l/E, 1}, {y, 


-1, 


1}, 


ContourStyle-> 


{Black} , 


PlotPoints->4 

] 

] 







Code Snippet 3. 




Figure A. 3: Output after running Code Snippet 2 then Code Snippet 3. 
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